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Abstract. Population studies of unidentified EGRET sources suggest that there exist at least three different 
populations of galactic gamma-ray sources. One of these populations is formed by young objects distributed along 
ON , the galactic plane with a strong concentration toward the inner spiral arms of the Galaxy. Variability, spectral and 

correlation analysis indicate that this population is not homogeneous. In particular, there is a subgroup of sources 
that display clear variability in their gamma-ray fluxes on timescales from days to months. Following the proposal 
by Kaufman Bernado et al. (2002), we suggest that this group of sources might be high-mass microquasars, i.e. 
accreting black holes or neutron stars with relativistic jets and early-type stellar companions. We present detailed 
inhomogeneous models for the gamma-ray emission of these systems that include both external and synchrotron 
. self-Compton interactions. We have included effects of interactions between the jet and all external photon fields 

to which it is exposed: companion star, accretion disk, and hot corona. We make broadband calculations to predict 
, the spectral energy distribution of the emission produced in the inner jet of these objects up to GeV energies. The 

results and predictions can be tested by present and future gamma-ray instruments like INTEGRAL, AGILE, 
and GLAST. 
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^ | 1. Introduction 120001 Bhattacharya et al. 12003(1 . These are bright sources 

r* ■ (isotropic luminosities in the range 10 34-36 erg s _1 ), with 

m > ■ Population studies of the unidentified gamma-ray sources an average photon spe ctral index T = 2.18±0.04 (F(E) oc 

KM I detected by the EGRET instrument of the Compton E -ry These SQurceSj whose number ig ^ 45±9j form thc 

b : Gam "! a : R 5 0bservator y ( 3rd EGRET catalog, Hartman Gamma _ Ray Population I (GRP I, Romero et al. EH®, 
et al. 11999(1 clearly suggest the existence of at least three 

different groups of galactic sources (Gehrels et al. 120001 Finally ' there is a third gr0Up ° f S0UrCeS that are dis " 

Grenier EBffl EES Romero EDS). There is a group of tributed forming a kind of halo around the galactic cen- 

weak and steady sources positionally correlated with the ter - These sources have hi S her luminosities, in the range 

Gould Belt (Grenier EMi Gehrels et al. TOU . These 1Q3 er § s • Thev have soft s P ectra ( r ~ 2 - 5 ) and dis " 

sources are thought to be nearby (100 - 300 pc), with P lav stron S variability (Nolan et al. 12303 ■ These sources 

isotro pic lum inosities - 10 33 erg s" 1 . Following Romero should be old ' with a S es from hundreds of Myr to Gyr. 

et al. (2001, we will call these sources the Local Gamma- About 45 ± 5 EGRET detections fall in this category 

Ray Population (LGRP). There are 45 ± 5 sources in this (Grenier E00D- They form the Gamma-Ray Population 

group. 11 ( GRP IT )- 

Another group of sources is concentrated along the Among GRP I sources there is a subgroup that displays 

Galactic Plane. They are well-correlated with star form- significant variability on time scales of weeks to months 

ing regions and HII regions, which is indicative of an as- (Torres et al.|2 001| N olan et al. 12003( 1 . Re cently , Kaufman 

sociation with young stellar objects (Romero et al. 119991 Bernado et al. $002) and Romero et al. I|2004(l have sug- 

Romero|200D- Log N - log S studies suggest that they are gested that this subgroup of GRP I sources might consist 

more abundant toward the inner spiral arms (Gehrels et al. of high-mass microquasars (i.e. microquasars formed by 

a compact object and an early- type stellar companion), 

Send offprint requests to: V. Bosch-Ramon where the gamma-ray emission arises from interactions be- 
e-mail: vbosch@am.ub.es tween relativistic particles in the jet and external photon 
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fields, most notably the stellar UV emission. Population 
studies by several authors (Grimm et al. 120021 Miyagi & 
Bhattacharva I2UIH|) suggest that high- mass X-ray bina- 
ries (HMXB) are distributed in the Galaxy following the 
spiral structure, presenting a steeper log AT - log S dis- 
tribution than low-mass X-ray binaries (LMXB), which 
means that LMXB are more uniformly distributed in the 
Galaxy than HMXB. Since high-mass microquasars form 
a subset of the whole set of HMXB, they are also expected 
to be distributed along the spiral arms, as is the case of 
GRP I sources. 

In the present paper we explore in more detail this hy- 
pothesis, presenting more realistic models for the gamma- 
ray emission. In particular, we will include effects of the 
interaction of the microquasar jet with the X-ray fields 
produced by the accretion disk and the hot corona that 
is thought to surround the compact object. We will also 
include synchrotron self-Compton emission, Klein-Nishina 
effects, and the back-reaction of the different losses in the 
particle spectrum of the jet. We will calculate, for some 
representative sets of the parameters that characterize 
high-mass microquasars, the broadband spectral energy 
distribution (SED) of these objects. 

We will start discussing the main phcnomcnological 
properties of GRP I sources in order to clarify what we 
should expect from our models. 

2. GRP I sources 

GRP I sources concentrate along the galactic plane and 
present a good spatial correlation with young stellar ob- 
jects (Romero et al. 1999). The variability analysis of these 
sources by Torres et al. I|2UU1|I clearly shows evidence for 
the existence of a subgroup with variable emission on 
timescales from weeks to months. This is corroborated by 
the recent results presented by Nolan et al. ( 2003), which 
are based on a maximum likelihood re-analysis of the 
EGRET data. These authors identify 17 variable sources 
within | 6° | from the galactic plane. These sources are 
clumped within | 55° | of the galactic center. 

A log A~ — log S analysis for all GRP I sources yields a 
distribution that is consistent with a power-law with in- 
dex (3 ~ 3.1 (Bhattacharya et al. 2003). This is far steeper 
than what is expected for a population uniformly dis- 
tributed along the galactic disk. For instance, for pulsars 
detected at 400 MHz the slope is (3 ~ 1.05. The unidenti- 
fied gamma-ray sources, on the contrary, seem to be con- 
centrated mainly in the inner spiral arms. To find possible 
further evidence for different populations among GRP I 
sources, we have implemented a log A" — logS analysis of 
both variable and non-variable low-latitude sources. 

First we have considered the 17 variable sources listed 
by Nolan et al. ( 20(21) • To take into account system- 
atic effects introduced by different exposure and back- 
ground resulting in non-uniform detectability, we have 
adopted the procedure described by Reimer (20TjT|l. The 
obtained log A" — log S plot is shown in Fig.^ lower panel. 
The normalized distribution can be fitted by a power- 
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Fig. 1. Log A" — logS 1 plot for gamma-ray sources within 
|6| < 6°. a) Upper panel: For non variable gamma-ray 
sources, b) Lower panel: For variable gamma-ray sources. 

law N(S) oc S- 13 , with [3 = 1.66 ± 0.31, significantly 
harder than for the entire sample. If we now consider those 
sources that classify as non- variable or dubious cases, we 
get the log N — log S plot shown also in Fig.f^ upper panel. 
In this case the distribution can be fitted by a power-law 
with index = 2.92 ± 0.36. 

The average spectral index is also different for both 
samples: in the case of the variable sources we have 
(r) = 2.04 ± 0.03, whereas for the remaining sources 
(r) = 2.16 ± 0.01. All this suggests that there are two 
different groups of sources, one formed by steady sources 
concentrated toward the inner spiral arms, and a second 
group with variable sources and a wider distribution along 
the galactic plane, although not as wide as that of radio 
pulsars (Bhattacharya et al. 12003(1 . 

Microquasars appear to be good candidates for com- 
pact and variable sources in the galactic plane. Since they 
can have very large proper motions (e.g. Ribo et al.[2002 
Mirabel et al. I2004|l . their distribution along the plane 
should be broader than that presented by supernova rem- 
nants and molecular clouds, which can be traced by star- 
forming regions and OB associations. Their spread, how- 
ever, is limited by the lifespan of the companion massive 
star, and hence it is not as extended as that of radio pul- 
sars. It is worth noting that Miyagi & Bhattacharya ( 2004 ) 
obtained for HMXB a (3 of about 1.9. In the next section 
we will discuss the potential of microquasars to be gamma- 
ray sources. 

3. Gamma-ray emission from microquasars 

Microquasars are variable sources at all wavelengths where 
they have been detected. At radio frequencies they present 
non-thermal jets and in some cases superluminal compo- 
nents (Mirabel & Rodriguez [1999). In the case of high- 
mass microquasars, with which we are concerned here, 
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the emission at optical and UV energies is dominated by 
the companion star, where it can reach luminosities of 
~ 10 39 erg s^ 1 . Accretion onto the compact object results 
in the formation of an accretion disk with typical tempera- 
tures of a few keV and thermal luminosities that can reach 
~ 10 37 erg s _1 . In the low- hard state, when the jets appear 
to be stronger, microquasars usually present a hard X-ray 
component that can be represented by a power-law plus 
an exponential cutoff at a few hundred keV. This compo- 
nent might be produced by a hot corona or ADAF region 
around the central object (e.g. Poutanen 1998, McClintock 
& Remillard 2004 ) although some authors have argued for 
a purely non-thermal origin in the jets (Markoff et al. 2001 
2003 . Georganopoulos et al. 120021) . Beyond the problem 
of the exact nature of the different contributions to the 
hard X-rays, it is clear that very relativistic electrons are 
present in the jets of microquasars. Since these electrons 
must traverse several photons fields, inverse Compton in- 
teractions are unavoidable. Such interactions can produce 
high-energy gamma-rays with a luminosity that will de- 
pend on the particular physical parameters that charac- 
terize the source. 

In recent years, several authors have studied mod- 
els for leptonic gamma-ray production in microquasars 
(Atoyan & Aharonian 119991 Paredes et al. 120001 Paredes 
et al. 2002, Kaufman Bernado et al. 2002, Romero et al. 
IM)2l Bosch- Ramon & Paredes E^tMll . Most of these mod- 
els consider steady and compact jets. Such jets are usually 
associated with the low-hard X-ray state of the sources. 
However, some microquasars present persistent radio jets 
and moderate X-ray emission. This is the case of the two 
already known microquasars that have been proposed to 
be counterparts of EGRET sources: LS 5039 (Paredes et 
al.EOnnillEI and LS I +61 303 (Massi et al-EDOl- The 
fact that a new microquasar candidate with similar X- 
ray properties has been recently discovered by Combi et 
al. pMjl within the EGRET location error box of 3EG 
J 1639-4702, lends additional support to the idea that per- 
haps high-energy gamma-ray emitting microquasars are 
not typical X-ray binary systems from the point of view 
of their X-ray behavior. 

In the present paper we will consider a microquasar 
with a continuous, persistent and inhomogeneous jet en- 
dowed with a structure similar to that adopted by Romero 
et al. (J2D03). In our case, however, we will focus on the 
leptonic content. The jet is ejected along the rotation axis 
of the compact object, which for simplicity is assumed to 
be perpendicular to the orbital plane 1 . 

The relativistic jet flow will move along the z axis with 
a bulk Lorentz factor Pj ct and a constant velocity [3 = v/c. 
Hadrons, although carrying most of the kinetic power, will 
play no radiative role in our model (for hadronic gamma- 
ray emission see the paper by Romero et al. 2003). We 
shall allow the jet to expand laterally, in such a way that 



1 Misaligned jets are probably a common case but our results 
will not be significantly affected by a small inclination angle. 
See Maccarone (2002). 



the radius R at a distance z from the compact object will 
be given by R{z) = £z 6 , with e < 1 and zq < z < z max . 
For e = 1 we have a conical jet. 

The relativistic electrons in the jet are assumed to have 
an energy distribution given by a power-law, as inferred 
from the observed synchrotron emission: 



iV(z,7e) = fc(z) 7c 



fc(*) = fc(z )(§ ) 2 = k(z )(^) 2e . 



CD 



(2) 



Here, p is the power-law index of the electron energy dis- 
tribution and Rq = R(zq) is the radius of the jet at the in- 
jection point. 7 C is the electron Lorentz factor and N(z, j e ) 
is the number of electrons of a given energy per unit of 
volume at a given distance from the compact object. The 
electron distribution does not depend on the direction of 
the particle velocity. 

Another important parameter associated with the elec- 
tron energy distribution is the maximum Lorentz factor, 
7emax- Its evolution is described by: 



7emax(z) — 7cmaxl+0j ^ 



(3) 



The parameter e is introduced to take into account loss or 
gain energy processes other than adiabatic losses, which 
are already counted through s, without a deep description 
of the involved physical processes (see, i.e., Ghisellini et al. 
1985). Both e and e cover our lack of knowledge of what 
happens to/within the jet. In the present work, we will set 
e to 1 from now on, to simplify the number of free param- 
eters, assuming as a first approximation a jet controlled 
by adiabatic evolution. We will come back to this issue in 
subsect. 15.11 analysing more carefully this question. 

Depending on the geometrical nature of the flow and 
assuming adiabatic expansion, the magnetic field changes 
with z as: 



B(z) = B(z Q ) 



= B 



<•)(?)• 



(4) 



Then, the synchrotron radiation density can be estimated 
as: 



Usyn,v{z) ~ T v (z) 



cC q 



\2dJ W^) 



-(p+4)/2 



(5) 



where s v is the source function of the synchrotron emission 
from an isotropic particle distribution and t„ is the syn- 
chrotron optical depth of the jet. The explicit expressions 
for both s„ and t v can be found in Pacholczyk (1970). 
We use the formulae in the optically thin case because 
the synchrotron emission mainly comes from frequencies 
beyond the self-absorption frequency, which is given by 



vx(z) m 2C 1 [R(z)C 6 k(z){m c c 2 )^-^}7T4B(z)^, 



(6) 
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Table 1. Constants in synchrotron formulae. 



Parameter (symbol) 


Value 




Ci 


6.27 x 10 18 eg 


s units 




1.37 x 1(T 23 c 


gs units 


C 6 (p = 2) 


8.61 x 10" 41 c 


gs units 



where Ci, C5 and Cq are numerical constants given in 
Tabled This is the local approximation to the synchrotron 
radiation field (Ghisellini et al. [1985). 

The total radiation field to which the leptons are ex- 
posed in the jet also will have a contribution from external 
sources. These contributions can be modeled as two black 
body components, one peaked at UV energies (the com- 
panion star field) and other at energies kT ~ 1 keV (the 
inner accretion disk field), plus a power-law with an ex- 
ponential cutoff at kT ^150 keV (the corona). With the 
exception of the disk, these contributions are assumed to 
be isotropic (Romero et al. 120021 Georganopoulos et al. 

EH. 

Once N(z, j c ) and the total energy density in seed pho- 
tons (C/ to tai(eo7 z)) have been determined, the IC interac- 
tion between them can be studied calculating the IC spec- 
tral energy distribution per energy unit. The cross section 
for both the Thomson and the Klein-Nishina regimes has 
been approximated by Blumcnthal & Gould (1970): 



da(x, £p,7c) _ 3o- T c 
de 4e 7| 

where 



/(*). 



(7) 



fix) = 



2x\nx + x + 1 - 2x 2 

x P(l/4 7c 2 , 1, x), 
c 



(4e 7 c x) 2 
2(1 + 4e 7cx) 



(1-x) 



4eo7e 2 (l-i)' 



(8) 
(9) 



eo and e are the energies of the incoming and the outgoing 
photons, respectively, ctt is the Thomson cross section, 
and 



P(l/4 7c 2 , 1, x) = 1, for l/4 7o 2 < x < 1, 



(10) 



and otherwise. 

The spectral energy distribution for the optically thin 
case in the jet's reference frame is: 



eL e — e 



:0max(z) flemaxf 



Omin(z) -'7emin( 



T,(z)U tot (e , z) 



xiV( 7c , z) d"/de a dz, 

de e 



(11) 



where S(z) is the surface of a perpendicular jet slice lo- 
cated at z. 

Notice that the external fields contributing to 
ftot(eo, z) should be transformed to the co-moving frame. 



Detailed expressions for such transformations are given by 
Dormer & Schlickeiser H2002H . In the observer's reference 
frame we have: 



e'L' = D 2+ r e 'L £ ,. 



(12) 



The integration is performed in the co-moving system and 
then the result is transformed to the observer's frame, 
hence the factor D 2+p , which is the Doppler boosting for 
a continuous jet. The energy of the scattered photon in 
the jet's reference frame (e) is boosted to e' = De. The 
Doppler factor D for the approaching jet is given by 



1 



r jet (1-/3 cos < 



(13) 



where (3 is the velocity of the jet in speed of light units and 
9 is the angle between the jet and the line of sight. We note 
that P(l/47 c 2 , 1, x) restricts the range of x to physical 
values, where e cannot be lower than e (x = I/47 2 ) or 
larger than e max (x = 1). 

In the case of the IC interactions with disk photons, 
a factor (1 — cos #)(p +1 )/ 2 must be introduced in Ea.f lllfl 
to take into account the fact that the photons come from 
behind the jet (Dermer et al.rj992). 

To make any calculation of the IC emission of a given 
microquasar, we have to specify first the jet power in lep- 
tons and hence the constant k in Eq. Q). In this work we 
shall adopt the disk/jet coupling hypothesis formulated 
by Falcke & Biermann l|1995l IT§9"9"|) . i.e. the jet power is 
proportional to the accretion rate: 



1$ = irR 2 cT jct f3 I m c c z k{zv) 1 L -Pd 1 



q c L a 



(14) 



Here, L a 



M acc c 2 is the accretion power onto the com- 



pact object and q c is a number < 1. The total jet power 
is Lj e t = qtot^acc = (<7hadrons + <?c)£ a cc- If the hadrons 
are also relativistic, ghadrons ~ <?tot >> 9c- Typically, 
<7tot ~ 10 _1 — 10 -3 and then q e is, in the kind of jet we 
are considering, in the range ~ 10~ 3 — 10~ 5 . 



4. Specific assumptions 

We make a number of specific assumptions regarding 
the characterization of the high-mass microquasar model 
adopted in the calculations of the SED. In particular, we 
will consider a system where the compact object is a black 
hole of 10 Mq that accretes 10 -8 Mq yr _1 . The gravita- 
tional radius R g for such an object is ~ 1.5 X 10 6 cm. The 
companion star has a radius of = 15 Rq and a bolo- 
metric luminosity of ~ 5 x 10 38 erg s _1 . Its radiation field 
is well represented by a black body with kT ^10 eV. The 
orbital radius is i? rb ~ 45 Rq. 

Around the black hole there is a thermal disk with 
typical temperatures given by kT ~ 1 keV. The total disk 
luminosity is ~ 10 37 erg s _1 . The hot corona above the 
central disk is represented by a power law with a hard 
photon index 1.6 and an exponential cutoff at ~ 150 keV. 
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Fig. 2. High-energy SED from IC up-scattering of stellar 
photons, a) Models with q c = 10~ 3 : for Fj C t =1.5 and 5. 
b) Models with q c = 10~ 5 : for r jc t =1.5 and 5. 



Fig. 3. High-energy SED from IC up-scattering of disk 
photons, a) Models with q c — 1CP 3 : for T jC t =1.5 and 5. 
b) Models with q c = 1(T 5 : for r jct =1.5 and 5. 



Regarding the jet, we assume it is injected at zq ~ 
50 R g . Its initial radius is Ro ~ 0.1 zq and its lateral ex- 
pansion is characterized by a coefficient e = 1. The elec- 
tron energy distribution has a power law index p = 2. The 
high-energy cutoff is assumed at 7 om ax(-Zo) = 10 , whereas 
the minimum Lorentz factor is 7 em in ~ 1- 

The magnetic field outside the coronal region, B(zq), 
is unknown. Most models for jet production require high 
magnetic fields, but these fields are usually attached to the 
inner accretion disk. Here we will assume a set of different 
possible values, in the range 0.1 — 200 G. These values are 
below equipartition, which could be reached further in the 
jet (e.g. in the radio emitting region). 

All these assumptions are summarized in Table El 

5. Calculations 

We have calculated the high-energy SED for a variety 
of cases using the typical values of the parameters given 
in Table |3 and different combinations of magnetic field 
B(zq), luminosity of the corona (£ cor ), bulk Lorentz fac- 
tor Tjot and disk/jet coupling parameter q c . In Table|31we 
give the range of these parameters adopted in our calcu- 
lations. 

In Figs. El to 03 we show some representative results for 
the SED produced by the IC up-scattering of stellar, disk, 
corona, and synchrotron photons, respectively. Curves ob- 
tained from different combinations of the parameters are 
shown in each panel. 

In Fig. El we show the whole SED of a representa- 
tive case, with the parameters listed in Table 0] We plot 
here the spectra for all contributions to the emission: the 
seed photon sources and the IC components. In Fig. we 
show a more extreme case, for a source with a high bulk 
Lorentz factor rj ct =10, a high-energy cutoff for electrons 
of 7emax(zn) = 10 6 , and a small viewing angle 9 = 1° (see 




6 7 
log (Photon energy [eV] ' 



Fig. 4. High-energy SED from IC up-scattering of corona 
photons, a) Models with q c = 10 -3 : for Lj Ct =1.5, 5; and 
L COI = 10 34 , 10 37 erg s" 1 . b) Models with q c = 10~ 5 : for 
L jct =1.5, 5; and L cnr = 10 34 , 



10 37 erg s 1 . 



Table ■ This would correspond to a microblazar, which 
is also a strong non-thermal X-ray source. It is difficult to 
say whether persistent microquasar jets could reach such 
high Lorentz factors, although it is an interesting case to 
explore. 

Finally, we present a 'realistic' case, reproduc- 
ing roughly the SED observed in the two EGRET 
sources that might be associated with known mi- 
croquasars: LS 5039/3EG J1824-1514 and LS I +61 
303/3EG J0241+6103 (see Fig. Eland Table EJ). For this 
particular case, we have taken both the disk and the 
corona to be faint, as it appears to be the case in both 
sources, with L cor = 3 x 10 32 erg s" 1 . Also, the mag- 



6 V. Bosch- Ramon et al.: High-mass microquasars and low- latitude gamma-ray sources 

Table 2. Fixed parameters in the models. 



Parameter (symbol) 


Value 


rilack hole mass (AZbh) 


i n A/T 
lUik/0 


Gravitational radius (-Hg) 


i.4o x iu cm 


Accretion luminosity (-Lace) 


in - 8 7i/f J2 „„„„ — 1 

1U Mqc year 


Stellar radius (7?*) 


10 JXq 


Stellar bolometric luminosity (£*) 


o x 1U erg s 


Viewing angle to jet's axis 6 


iu 


Distance from jet's apex to the compact object (-^o) 


OU xtg 


Initial jet radius (Rq) 


U.izo 


Orbital radius (i? rt>) 


40 it© 


Peak energy of the disk ffeT^jsk) 


i Kev 


Peak energy of the corona 


loU kev 


Peak energy of the star (fcT star ) 


10 eV 


Expansion coefficient of the jet (e) 


1 


Minimum Lorentz factor for electrons in the jet (jet frame) (7emin) 


~ 1 


Maximum Lorentz factor for electrons in the jet (jet frame) (7emax) 


10 4 


Electron energy distribution power-law index (p) 


2 


Photon index for the corona (r c0 r) 


1.6 


Total disk luminosity (Ldisk) 


10 37 erg s" 1 
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Fig. 5. High-energy SED from IC up-scattering of syn- 
chrotron photons, a) Models with q c = 10 -3 : for Tj ct =1.5, 
5; and B(z )=0.l, 100 G. b) Models with q c = 10~ 5 : for 
r jct =1.5, 5; and B(z o )=0.l, 100 G. 



Table 3. Range of parameters adopted for calculation of 
different models. 



Parameter (symbol) 


Value 


Corona luminosity (L COI ) 


10 34 , 10 37 erg s" 1 


Magnetic field (B(z )) 


0.1, 100 G 


Bulk Lorentz factor of the jet (Fj e t) 


1.5, 5 


Disk/jet coupling constant (q B ) 


10~ 5 , 10 -3 



Table 4. Parameters for a representative microquasar 
model. 



Parameter (symbol) Value 



Corona luminosity (L C or) 10 35 erg s 1 
Magnetic field (B(z )) 10 G 

Maximum electron Lorentz factor (-y C m&x(zo)) 10 4 
Bulk Lorentz factor of the jet (lj e t) 2.5 
Disk/jet coupling constant (g c ) 10 -4 



netic field required now to match the observations is 
i~ 200 G, q e is ~ 10 -3 , and consistent with observational 
constraints from both sources, a mildly relativistic jet ve- 
locity (rjct = 1-1) is used. The electron high-energy cutoff 
is at 7omax(-2 ; o) ~ 10 4 , in agreement with EGRET data 
(see Fig. El . 

5.1. Electron energy evolution 

In this subsection, we discuss whether our parameteriza- 
tion is consistent with the local cooling rates of the elec- 



trons due to energy losses. This will allow us to know to 
what extent the evolution of the electrons is well described 
through the adopted parameterization, that corresponds 
to adiabatic expansion. If radiative losses are smaller than 
adiabatic ones, then the parameterization is fine at least 
as a first order of approximation. 

We have determined quantitatively the importance of 
the energy density of the different seed photon fields, tak- 
ing into account in our computations both the internal 
ones (synchrotron photon and magnetic fields) and the 
external ones (star, disk and corona photon fields). In 
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Table 5. Parameters for an extreme case (microblazar). 



Parameter (symbol) 



Value 



Corona luminosity (L COI ) 
Magnetic field (B(z )) 
Maximum electron Lorentz factor (7 e 
Jet bulk Lorentz factor (rj ct ) 
Viewing angle to jet's axis 9 
Disk/jet coupling constant (q e ) 



10 



<(*>)) 



35 erg s 1 
10 G 
10 6 
10 
1° 

io~ 4 



Table 6. Parameters for a 'realistic' case. 



Parameter (symbol) 




Value 



2 4 6 

log (Photon energy [eV] ) 



Corona luminosity (L COI ) 
Magnetic field (B(z )) 
Maximum electron Lorentz factor (7 e; 
Jet bulk Lorentz factor (Pj e t) 
Viewing angle to jet's axis 6 
Disk/jet coupling constant (q e ) 



3xl0 32 erg s" 
200 G 
10 4 
1.1 
10° 
10" 3 




-1 1 2345678910 
log (Photon energy [eV] ) 



Fig. 6. Complete SED for a model with q e 



1CT 



B(z Q )=10 G, L COI = 10 35 erg s"\ 7cmas (z fl ) = 10 4 , 
rj et =2.5, and a viewing angle of 10°. This is the repre- 
sentative case. 



Table [7J we show the energy densities for the mentioned 
photon fields in the comoving frame at the base of the jet 
for the three microquasar models we have calculated (see 
Tables 0113 and EJ. In Fig. 03 we show the evolution of the 
cooling times along the jet for the adiabatic and the radia- 
tive losses. This allows us to find out which one dominates 
at different distances. To compute and compare the differ- 
ent cooling times, we have taken the electron Lorentz fac- 
tor which would correspond in our model to a given value 



Fig. 7. SED for a model with q c = 1CT 4 , B(z )=10 G, 
icor = 10 35 erg s~\ 7cmax(^o) = 10 6 , r jot =10, and a 
viewing angle of 1°. This is the extreme case. 



IC emission 

seed photons 

EGRET sources 




-4 -3 -2 -1 1 2 3 4 5 6 7 8 9 
log (Photon energy [eV] ) 



Fig. 8. SED for a model with q c = 1CT 3 , B(z )=200 G, 
Lcor = 3 x 10 32 erg s -1 , 7cmax(zn) = 10 4 , r je t=l-l, and a 
viewing angle of 10°. This is the 'realistic' case. 



of z, considering the adiabatic evolution of electrons. In 
this plot, we show the behavior of the 'realistic' and the 
representative cases as examples; the microblazar case is 
not too different. 

From Fig. it is clear that it is enough in the 'real- 
istic' case to account only for adiabatic losses to describe 
the electron energy evolution along the jet, whereas it is 
not sufficient for the representative case at small values 
of z. In the second situation, however, this is not a prob- 
lem since we do not pretend to describe the microphysics 
within the jet, but only to parametrize it in such a way 
that the expected emission will be observed. This parame- 
terization could take into account reacceleration processes 
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Table 7. Energy densities at the base of the jet in units of erg cm~ 3 . 



Model 


Jet velocity 


Lorentz factor 


magnetic 


synchrotron 


star 


disk 


corona 




(c) 




field 


photon field 


photon field 


photon field 


photon field 


Representative 


0.917c 


2.5 


2 


3.2 


220 


2.4xl0 5 


2.2xl0 6 


Microblazar 


0.995c 


10 


2 


49.8 


1.6xl0 4 


3.7 


3.4xl0 4 


'Realistic' 


0.417c 


1.1 


800 


6.3xl0 4 


7.6 


87.8 


8.2xl0 4 



10= 




a) 


— adiabatic losses 

— radiative losses 




star losses 


10 3 






corona losses 






10 1 












1Cf 1 












10 5 




l» 








10 3 






— adiabatic losses 

— radiative losses 




star losses 


10 1 






corona losses ^ 






10"' 












io- 3 












10° 




10 9 


10 


10 



z [cm] 

Fig. 9. Cooling time evolution with z for adiabatic (solid 
line) and radiative (long-dashed line) losses. Two re- 
gions in the plot have been established depending on the 
dominant source of radiative losses, which happen to be 
the corona IC losses close to the compact object and 
the star IC losses farther away, a) 'Realistic' case, b) 
Representative case. 

acting on electrons. The assumption that electrons are ac- 
celerated along the jet is natural for many jet models (see, 
i.e., Biermann & Strittmater 1987). Future work will be 
done to implement this in a more explicit way. 

6. Comments 

We can see from Figs. 12151 that for several models with 
disk/jet coupling constant q c = 10~ 3 we can get the ex- 
pected luminosities in the observer's frame inferred for 
GRP I sources with the right photon index at energies 
~ 1 GeV, i.e. r ~ 2. When the magnetic field is strong 
enough (B(zq) ~ 100 G), SSC emission alone can account 
for luminosities ~ 10 35 erg s _1 at 1 GeV. Models with 
bulk Lorentz factors Tj Ct =1.5 and 5 do not produce dra- 
matically different results for a jet with a viewing angle 
of ~ 10° except for the case of the scattered stellar pho- 
tons. In case of lower magnetic fields, the IC scattering 
upon external fields clearly dominates. For instance, when 
B(zq) = 0.1 G, the IC component upon stellar seed pho- 



tons is more than 5 orders of magnitude greater than the 
SSC emission at 100 MeV. For a similar value of q e and 
a Fjet = 5, the stellar IC component usually dominates 
over the disk and corona components. Only in the case of 
a powerful corona contribution (L COI ~ 10 37 erg s _1 ) can 
the latter surpass the IC emission produced in the stellar 
field. The contribution arising from interactions with disk 
photons reaches values ~ 10 34 erg s _1 only for energetic 
jets (q c = IO" 3 ). 

Models with light leptonic jets (q c = 10 -5 ) can pro- 
duce significant gamma-ray sources (~ 10 34 erg s _1 at 
100 MeV) when a strong corona is present (see Fig. 0J. 
In these models, the spectrum tends to be a bit harder 
than in the case of stellar photons, with our current set of 
assumptions (value of r cor , etc). 

Comparing Figs. and it is possible to see the dif- 
ferences between a 'mild' microquasar, with Fj et = 2.5, 
a viewing angle of 9 = 10° and q = 10~ 4 , and a more 
'extreme' case with Lj t — 10 and 9 — 1°, i.e. a microb- 
lazar. Since we have assumed in both cases a magnetic 
field of B(zo) = 10 G, external IC scattering always domi- 
nates over SSC. In the 'mild' microquasar the emission at 
MeV-GeV energies is mainly due to upscattering of coro- 
nal photons whereas, in the extreme case, the IC emission 
in the stellar field exceeds by far the other contributions. 
The hard X-ray counterpart can be important, even be- 
yond the cutoff for the corona (at ~ 150 keV). Hence, 
INTEGRAL observations can be useful to unmask candi- 
dates that are obscured at optical wavelength (see Combi 
et al. 120041 for a recent study in this direction) . 

An interesting feature of the extreme case is that the 
gamma-ray emission with a hard spectrum P ~ 2 ex- 
tends up to tens of GeV. At TeV energies, the spectrum 
becomes very soft due to the Klein-Nishina effect. This 
type of sources should be detectable with modern Imaging 
Cherenkov Telescopes like HESS and MAGIC. Even sys- 
tems with low-mass companions and weak coronas might 
be detectable since the SSC component has luminosities 
of ~ 10 35 erg s _1 . 

In the case we have called 'realistic' (Fig.|HJ), the high- 
energy emission is dominated by SSC. We emphasize that 
this case is more realistic than the previous ones only 
in the sense that it fits well the broadband spectra of 
known microquasars that are suspected to be high-energy 
gamma-ray sources. A basic feature of these objects is 
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their weak X-ray luminosity, a feature that is shared by 
most of the unidentified EGRET sources. 

We now comment on the main differences between the 
models presented here and those already published by 
Georganopoulos et al. ( 2002 ) and Bosch- Ramon & Paredes 
(2004). The first authors consider just the interaction of a 
homogeneous jet with disk and stellar photon fields. They 
work in the Thomson regime with the aim of producing the 
hard X-ray emission observed in objects similar to Cygnus 
X-l as a result of the IC contributions. No magnetic field is 
considered and adiabatic losses are not taken into account. 
On the contrary, we adopt a thermal Comptonization field 
at hard X-rays, possibly originating in a corona around the 
black hole (e.g. Zdziarski et al. l2UU3f) . as suggested by the 
observation of a Compton reflection feature and an iron 
Ka line in some sources (for an alternative interpretation, 
see Markoff et al. 200 3b"|) . In this sense, our treatment is 
more similar to that of Romero et al. lP002). However, we 
include the effects of the magnetic field as well as a more 
sophisticated parameterization of the jet. SSC emission, in 
fact, seems to dominate for some reasonable choices of the 
parameters. In addition, wc go to much higher energies, 
taking into account the Klein- Nishina effect. Regarding 
the Bosch-Ramon & Paredes f2QDlt model, here we intro- 
duce several improvements, from the corona effect up to a 
z-dependence of the magnetic field and other parameters. 

7. Discussion 

Variability in the gamma-ray emission of microquasars on 
timescales from days to months can be caused by both 
changes in factors external to the jet or in the jet itself. 
The winds of the companion star can change inducing vari- 
ations in the accretion rate (Reig et al. 120031 McSwain 
et al. I2004J) . Also, if the orbit of the compact object is 
not completely circularized, periodic changes in the IC 
flux from the stellar photon field and in the accretion rate 
can be expected. These changes would vary, in turn, the 
power carried by the jet (through changes in q c ), the emis- 
sion of the jet, and the emission of the thermal plasma 
around the black hole. Precession of the accretion disk 
(with timescales of months, e.g. Brocksopp et al. I1999JI 
might also result in a time-modulation of the correspond- 
ing soft X-ray field. Precession of the jet itself can, ad- 
ditionally, produce strong variability due to the variation 
of the Doppler factor (Kaufman Bcrnado et al. l2002J) . All 
these contributions might result in a complex lightcurve, 
with several associated timescales. 

In sources where the jet is steady and a transient fea- 
ture associated with the low-hard state, its disappearance 
in the high-soft state would lead, of course, to the sup- 
pression of the gamma-ray emission. This should be an 
additional source of variability. This transition, however, 
seems not to occur in microquasars with persistent jets 
like LS 5039. 

Besides medium and long-term variations, the presence 
of shocks in the jets can introduce very rapid changes in 
the gamma-ray flux. Relativistic shocks are the natural 



result of sudden changes in the injection rate of plasma in 
the jets. These shocks propagate downstream increasing 
the energy of the particles in the fluid and amplifying the 
magnetic field. When the shock finds a small feature in 
the underlying jet (e.g. an inhomogeneity in the particle 
density or a bend in the flow direction, see Romero|1996) 
a very rapid variation in the high-energy flux can occur. 
The typical timescale for these events will be determined 
by the time it takes the shock to move trough the feature. 
Assuming that the feature has a size of the order of the 
jet radius, this gives: 

^var ~ "pr j (15) 

where D s is the Doppler factor of the shock and v s its 
velocity. For i?j et ~ 10 7 cm, and a relativistic shock with 
D s »j 5 and v s ~ c, we have t val ~ I0 -4 s, so this would be 
a very rapid flickering superposed on the longer variability. 

Microquasars are not, of course, the only kind of 
galactic object that might display variable gamma-ray 
emission. Other alternatives include early-type bina- 
ries (Benaglia & Romero 2003), accreting neutron stars 
(Romero et al. I2()01|l , pulsar wind nebulae (Roberts et al. 
2002) and exotic objects (Punsly et al. 2000). However, mi- 
croquasars are perhaps the most attractive candidates to 
explain a significant fraction of the variable GRP I sources 
because of the presence of relativistic jets in these ob- 
jects, as well as the external photon fields provided by the 
companion star and the accreting matter. GRP II sources 
might also be associated with old low-mass microquasars if 
the jets of these objects have magnetic fields strong enough 
to make SSC the primary source of radiation (Kaufman 
Bernado et al. l2"UDl)) . 

8. Conclusion 

We have shown that the variable gamma-ray sources found 
on the galactic plane have some common features that 
make it reasonable to consider them as a distinctive group 
of GRP I sources. We have suggested that these sources 
might be microquasars with high-mass stellar compan- 
ions and we have developed some detailed models to ex- 
plain the gamma-ray production in this type of objects. 
In particular, we have considered inhomogeneous jet mod- 
els where gamma-rays are the result of inverse Compton 
interactions of leptons in the jet with locally produced 
synchrotron photons as well as external photon fields. We 
have calculated the emission resulting from the upscatter- 
ing of disk, coronal, and stellar photons, incorporating a 
full Klein-Nishina calculation and the effect of losses in 
the particle spectrum. We have shown that a wide variety 
of spectral energy distributions at high energies can be ob- 
tained from different and reasonable combinations of the 
physical parameters like magnetic field, jet power, coronal 
and disk luminosities, etc. It seems clear that the micro- 
quasar phenomenon can be naturally extended up to the 
highest energies and that we can expect these objects to 
manifest themselves as a distinctive group of gamma-ray 
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sources that might be detectable with satellite-borne in- 
struments like those to be carried by AGILE and GLAST, 
and even by ground-based Cherenkov telescopes like HESS 
and MAGIC. 
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